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Abstract 

We present a preliminary measurement of Ri, R2, and , which are the three parameters used 
to characterize the B^D*iU(, form factors (^1, V , and A2). We use 85 x 10^ S^-pairs accumulated 
on the T(4S) resonance at PEP-II. In this analysis we use the decay mode B^^D^*e~ v and its 
charge conjugate. The Z?*"*" is reconstructed in the channel D*^ —>■ D^tt^ and the in the 
channel — > K~Tr~^. We parameterize the form factors in terms of their ratios (determined by 
the parameters Ri and R2) and the common slope in the variable w (a quantity related to 
the momentum transfer in the decay process). These parameters are determined via an unbinned 
maximum likelihood fit to the distributions in four kinematic variables (three decay angles and w). 
The resuhs are Ri = 1.328 ± 0.055 ± 0.025 ± 0.025 and R2 = 0.920 ± 0.044 ± 0.020 ± 0.013 for the 
ratios and = 0.769 ± 0.039 ± 0.019 ± 0.032 for the slope. The stated errors are the statistical 
uncertainty from the data, statistical uncertainty from the size of the Monte Carlo sample and the 
systematic uncertainty, respectively. 
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1 Introduction 



The decay B^D*lv can be described by three form factors: two axial form factors, Ai and A2, 
and one vector form factor, V . They are functions of the momentum transfer (f' (or equivalently, 
u;, defined below). Measurements of these form factors provide us with a "laboratory" for testing 
the predictions of heavy quark effective theory (HQET) ^ . These form factors are related to each 
other by heavy quark symmetry (HQS) through the formalism of heavy quark effective theory. 
Deviations from the HQET relationships can be computed as corrections to the theory. They can 
also in principle be measured; the parameters we adopt for this analysis are inspired by HQET, 
but allow for such deviations. 

We introduce here a novel method of extracting these parameters from the data. We use 
an unbinned-maximum-likelihood method, but introduce approximations that allow us to correct 
efficiency and resolution with the limited Monte Carlo (MC) data sample available. In Sec. El we 
sketch out how these approximations work and how we evaluate their impact on the errors. 

Improved measurements of the form factor parameters R2 and yield a significant reduc- 
tion in the systematic error obtainable on the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
Vcb- In this analysis, though we restrict ourselves to the cleanest available decay channel, we still 
obtain nearly 17 times as many reconstructed B^D*lu decay candidates as the pioneering CLEO 
analysis P]. This already produces a substantial improvement in the error achieved on the param- 
eters Ri and R2 needed for extracting Vcb- Further, a better understanding of this channel, which 
is the dominant background to i? — > Xuiv, is also of great importance to the study of inclusive and 
exclusive B Xuii^ decays and the determination of Vub- 

With the addition of more data and the use of the other available decay modes we will be able 
to further probe the consistency of HQET predictions in the near future. 

2 Formalism 

The Feynman diagram for the decay B^ — > D*~^i~V is shown in Fig. ^ 



The square of momentum transfer from the B to the D*, g^, is linearly related to another 
Lorentz invariant {w) by 




i 



d D 



Figure 1: Quark level diagram that leads to the decay B D*lu 



w = 




= VB ■ VD* = 



PB ■ PD* 

MbMd' 



(1) 



2MbMd* 
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where Mb and Mjy* are the masses of the B and the D* mesons, ps and pd* are their four 
momenta, and vb and vd* are their four- velocities. In the B rest frame w reduces to the Lorentz 
boost factor 7/)* = (Eu* /Md*). 

The range of w and are restricted by the kinematics of the decay with = corresponding 

to 

_ M| + Ml. 



2MbMd* 



1.504 



(2) 



and Wr, 



1 corresponding to 



'imax 



,2^2 



(Ms - Md'Y - 10.69(GeV/c^) 



(3) 



2.1 Kinematic Variables 

We choose to reconstruct the mode where the D* decays to Dtt, and the -D to a Kit . The outgoing 
particles are shown in Fig. |21 The pion from the D* decay and the electron are directly detected, 
and the D is reconstructed from its daughters, the K and vr. 




Figure 2: Kinematics of B^D*lv decays. 



In Fig. IJJwe define three angles, which, along with w, constitute the four independent kinematic 
variables we use to characterize this decay: 

• cos 9 1, the cosine of the angle between the direction of the lepton in the virtual W rest frame, 
and the direction of the virtual W in the B rest frame 

• cos the cosine of the angle between the direction of the D in the D* rest frame, and the 
direction of the D* in the B rest frame 

• Xi the azimuthal angle between the plane formed by the D* — D system and the plane formed 
by the W — i system. 

A D*lu decay is completely characterized by the four kinematic variables w, cos Of, cosOy and 
X- Since the B meson is spinless, the direction of the D* — W axis relative to the B direction is 
irrelevant to the dynamics. 
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2.2 Four Dimensional Decay Distribution 



2.2.1 Helicity amplitudes 

The Lorentz structure of the B^D*!}/ decay amplitude can be expressed in terms of three am- 
plitudes -ff-, and Hq)^ which correspond to the three allowed polarization states of the D* 
(two transverse and one longitudinal). These amplitudes can be completely specified in terms of 
the axial and vector form factors as follows pQ: 



H+{w) . 

1 



-{Mb + Md*)A^{w) + 2-^^^^^V{w), 



-{Mb + Md*)Ai{w) - 2 



Mb + Md- 
Pd*Mb 



2Md' 



Mb + Md 
Ai{w){Mb + Md'){MI - Ml, - g2 



V{w), 



(4) 
(5) 
(6) 



Mb + Md* 



A2{w)) 



where pjj = Md* ^/w"^ — 1 is the magnitude of the momentum of the D* in the B rest frame. 

By contracting the relevant lepton and hadron tensors, neglecting small terms of order m^/q^, 
and doing the phase space integrations, the full differential decay rate 



dT{B^D*iv) _ &Gl\Vcb\'^MBMl,ryJ{w^ - - 2-wr + r^) 

X 



dw dcos 9i dcos 9v dx 8(47r)^ 

H+^{l-coseif sin^ 9v + HJ{l + coseif sin^ By 
+4Ho'^ sin^ Oi cos^ By - 2H+H_ sin^ Oe sin^ Oy cos(2x) 
—4:H+HQsm6i{l — cos 0f)sin 0ycos 0y cos x (7) 
+4:H^HQsm6i{l + cos 0^ )sin Oycos 9y cos x] 

can be obtained in terms of the helicity amplitudes. Here r = (y ^lf* ^ is the ratio of the D* to the 
B mass. Details of the derivation can be found in Neubert^. The four-dimensional distribution 
of w, cos 9 1 cos 9y and x that is described by Eq.ljT)) is the physical observable from which we can 
extract the form factors. 

2.2.2 Heavy Quark Symmetry Relationships 

HQS relates the three form factors to each other, as follows: 



R2{w) 2 

i?*2 w + 1 
Ri{w) 2 

^ 1 

where we have defined the constant 



^2H = ^77^^! (^), (8) 
nw) = '-^S^^A,{w), (9) 
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Ri{w) = R2i'w) = 1.0, i.e., perfect HQS, implies that the form factors, A2 and V are identical for all 
values of w. Since HQS is not exact, Ri and R2 can differ from 1.0 and exhibit some t/;-dependence. 
Ai can be related to the Isgur-Wise function by 

A,{w)=R*^^hAAw), (11) 

where in the HQET limit is the Isgur-Wise function £,{w) 0. HQET predicts /i^^(l) = ^(1) = 
1.0. 

It is convenient to re-express the helicity amplitudes Hi [i = ±,0) in terms of the functions 
Ri{w) and R2{w) as follows 



Hr 



H. 



v/(l - 2wr + r^) 1 



w — \ 
w + l 



Ri{w] 



1 + 



(l-r) 

{w-l){l-R2{w)) 
1 — r 



(12) 



V y w + 1 

(l-r) 



where the helicity amplitudes that appear in Eqs. (|4l6j) differ from these only by a common kinematic 
factor times /iaj : 



-Mb 



R*{l-r'^){w + l) 
2VI - 2wr + r'^ 



(13) 



The function (w) can be conveniently parameterized by exploiting the fact that w — 1 is a small 
parameter ( < 0.5) to express it as a power series expansion 



hA.iw) = hA, (1) (1 -p^iw-l)+ c{w - if + ...) 



(14) 



where p'^ is called the slope and c is called the curvature (note that in the literature p'^ is often 
referred to as /O^J. Corrections to HQET modify /lAi(l) and thus lead to deviations from the 
HQET limit of Hai = 1.0. However, in this analysis we only deal with the shape and relative 
normalization of the form factors, and consequently the overall normalization is irrelevant. 

For this preliminary analysis we only use the expansion to first order in {w — 1), so we set c = 0. 
In our baseline analysis we treat Ri and R2 as constant. However, we also show results in which 
we extend the analysis to explore the impact of w-dependence. 



2.3 Form Factor Predictions 

As discussed in Sec. 12.2.21 for infinitely massive b and c quarks, we expect Ri = R2 = 1.0, but for 
finite masses they are modified by both perturbative and non-perturbative corrections. 
Calculating higher-order loop corrections to the form factors yields expansions of the form: 
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R^{w) = 1 + [a, (...) + a2 (...)] + (^) (...), (15) 

R2{w) = 1 + [a, (...)+ a2 (...)] + (^) (...). (16) 

The coefficients of the as terms (shown as ehipses) have been calculated perturbatively up 
to second order, which gives confidence that they are accurate to about one percent (see 0). 
The coefficients of the terms are combinations of quantities called "subleading Isgur-Wise 

functions." 

Different models in the HQET framework evaluate the subleading Isgur-Wise function correction 
terms differently, resulting in a variety of predictions for Ri{w) and R2{w). Neubert based on 
work with collaborators, in the early 90's predicted 



Riiw) = 1.35 - 0.22(u' - 1) +0.09(?i; - l)^ (17) 
R2iw) = 0.79 + 0.15(t(; - 1) -0.04(w; - 1)2. (18) 

More recently, Caprini et.al. using spectral functions, dispersion relations, and heavy quark 
symmetry to evaluate the non-perturbative terms predict 



Riiw) = 1.27- 0.12(u' - 1) +0.05(w - l)^ (19) 
i?2(u>) = 0m + 0.11{w-l)-0.06{w-lf. (20) 

Ligeti and Grinstein [5] using similar HQET machinery find 



Ri{w) = 1.25 - 0.10(u; - 1), (21) 
R2{w) = 0.81 + Om{w - 1). (22) 

Whereas all the above HQET-based predictions predict form factor values within a fairly narrow 
range, older predictions relying on ad-hoc potential models vary widely, e.g., Close & Wambach 
using a simple quark model predict [HI 



Ri{w) = 1.15 - 0.07(u; - 1), (23) 
R2{w) = 0.91 + 0.04(u; - 1). (24) 

It can be seen that in all the above predictions the coefficients of the {w — 1) and {w — 1)^ terms 
are quite small because Ri and R2 are by construction as ratios expected to vary only slightly 
with w, where /i^i has no such restriction. 



3 The BaBAR Detector and Dataset 

The data set used in this analysis was collected with the BABAR detector at the PEP-II storage 
ring during the period between 2000 and 2002. It corresponds to 84 fb~^ collected on the T(4S) 
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resonance, which yields approximately 85 x 10^ BB-pairs. There are ~ 8.5 x 10^ B^ D*lv 
decays in this sample of which we have reconstructed 16,386 B^^D*^e~h'e (or charge conjugate) 
candidates using only the D^tt^ decay of the D* and the K~it~^ decay mode of the . 

The BABMi detector is described elsewhere in detail This analysis uses four of the five 
subdetectors oi BABMi: the silicon vertex tracker, the drift chamber, a Cerenkov-light-based particle 
identification detector, and the electromagnetic calorimeter. This analysis depends critically on the 
silicon vertex tracker to reconstruct the low momentum pions produced by the decay D*^ D^tt^, 
about two-thirds of which do not traverse more than about a fourth of the drift chamber. 

4 Reconstruction and Event Selection 

We reconstruct the lepton, and the D* through its decay products. After reconstruction of all 
tracks, we choose cuts that select B^D*iv decays and determine the kinematic variables w, cosOg, 
cos 9v and x- 

4.1 Event selection 

In this analysis we only consider the decay channel D*~^e~i'e with the D*~^ decaying to D^tt~^ and 
the decaying in the K^ir'^ mode. In case of multiple D*ev candidates in a given event we 
choose the candidate with Ktt mass {mKir) closest to the mass. 

For event selection we use the procedure developed for our Vcb analysis [S], and we also apply 
the same selection criteria. We summarize the most salient of these cuts here: 

• The momentum of the lepton in the center-of-mass (CM.) frame is required to be larger 
than 1.2GeV/c. This criterion selects B semi-leptonic decays and suppresses continuum (cc) 
and cascade [b^c^t) backgrounds. 

• The momentum of the D* in the CM. frame must be between 0.5 and 2.5 GeV/c. 

• The soft pion from the D* decay must have a transverse momentum greater than 50 MeV/ c. 

• The probability of the fit of the D*l vertex to the beamspot constraint must be greater 
than 1%. 

• To further suppress continuum background, we select only candidates with | cos ^thrust I < 0.85, 
where ^thrust is the angle between the thrust axis of the D*lv candidate and the thrust axis 
of the rest of the event. 

The cosine of angle 9 by between the direction of the B and the direction of the D* — I system 
can be computed from the kinematics of the decay (see next Section). Candidates with cosOby 
between —10 and +6 have been used to estimate background. We include only events that have 
I cos I ^ 1-2 in the final sample. 

The final selection is based on 6m = mKnns ~ n^'Kn- For candidates in which the vr^ track is 
reconstructed with twelve or more drift chamber hits we cut on the mass difference 0.144 < 5m < 
0.147 GeV/c. For the less well-measured case with fewer than twelve drift chamber hits we loosen 
the cut to 0.143 < 6m < 0.148 GeV/c. About three-quarters of the candidates are in the latter 
category. 
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4.2 Kinematic variable determination 



Without the neutrino we do not have sufficient information to fully reconstruct the kinematic 
variables w, cos Oi, cos By and x- However, we can use the measured variables to construct an 
adequate approximation. Using energy-momentum conservation and assuming that the neutrino 
mass nil, is zero, we have 

= ml = Ml + M^- 2EbEy + 2pBpY cos 6 by, (25) 

where pY = Pf) + Pe is the four momentum of the combined D* and lepton, My = py is the mass 
squared and py is the three-momentum. The B meson energy Eb and three-momentum pB can be 
estimated from the energies of the colliding beam particles, so we can solve for cos 9 by as follows 

M| + M^ - 2EbEy 
cos e BY = ^-^ ■ (26) 

Thus we can determine the angle between the B and the direction {Y = py/py) of the D* — i 
system, but we do not know its azimuthal angle (f>BY around this direction. This is illustrated in 
Fig. 13 The direction of the B must lie on the cone around Y defined by the angle Oby- 

For each possible 4>by we can compute the kinematic variables w, cos 9£, cos By, x- Since we do 
not know which direction is correct we perform an average to estimate each variable. We average 
over four points: two in the D* -lepton plane corresponding to the azimuthal angles (pBY = and vr 
and two points out of the plane corresponding to the angles 12. Further, since BB production 
follows a sin^ 9b distribution in the angle between the B direction and the beam collision axis, we 
weight the kinematic variables evaluated at each point by the value sin^ 9 b corresponding to the B 
direction for that point. 

Fig.|l]illustrates the resolution achieved by this 'partial reconstruction' technique. It is seen that 
the core widths for each resolution distribution are relatively small compared to the full width of 
each kinematic variable. The resolution is dominated by the average over the B direction; detector 
resolution makes only a relatively minor contribution. The low-side tail on cos 9ii can be attributed 
to final state radiation. 

The resolutions of the four kinematic variables are highly correlated. A simple factorized reso- 
lution function in which the resolution is represented by a product of independent function for each 
variable fails to capture these correrlations. Thus, we are dependent on the Monte Carlo simulation 
to account for resolution effects. 

The distributions of the reconstructed kinematic variables cos 9i, cos9v and x are displayed 
in Fig. [3 The shaded region is the distribution of the background as estimated from the MC 
simulation using the method described in Sec. IHl below. It can be observed that the background 
distribution is much smaller than the signal contribution. 

5 Simulation 

This analysis is dependent on Monte Carlo (MC) simulation to model the efficiency and the back- 
ground distributions. The degree to which we can trust our simulation to model both the detector 
and the underlying physics processes largely determines our systematic errors. 

For the detector simulation we use BABAR's GEANT4 based simulation jHI- The simulation has 
been validated and sophisticated correction procedures have been developed using many control 
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Figure 4: MC assessment of the experimental resolution for the variables cosBi, cosOy and x- 
For each variable we plot the difference beween reconstructed and generated values. The resolution 
is generally small compared to the ranges of the variables as shown in Fig. |SJ 
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Figure 5: Distribution of kinematic variables w, cosOg, cos By and x foi^ selected 
events with estimated backgrounds (shaded region). 

samples. For simplicity, we do not apply these correction factors, but use them to evaluate sys- 
tematic errors (which turn out to be small for this analysis, thereby justifying using them for a 
posteriori error estimation rather than as correction factors). Event generation and particle decay 
are modeled using the package EvtGen |lUj . 

5.1 Slow pions 

Of particular importance to this analysis is the modeling of the efficiency for detecting low- 
momentum pions. This is a difficult task since low momentum pions are lost through the interplay 
of acceptance, decay-in-flight, and stopping and scattering in the beam pipe or vertex detector. 
The performance of our simulation can be seen in Fig. ^ which shows the distribution of cosine 
of the helicity angle (cos Oy) for D*~^ — > D^n^ obtained from inclusively produced D* mesons for 
both the data and the MC sample in bins of center-of-mass momentum. The MC sample has been 
corrected by a factor of the form 

fcorr = + a cos ev + 6 cos^0v)) (27) 
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Figure 6: Distributions of the cosine of the hehcity angle cos Oy for data (points) and MC (his- 
togram) in bins of D* center-of-mass momentum. The momentum range for each plot is given its 
upper left-hand corner. The MC has been fit to the data as described in the text. 
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where the normahzation and the parameters a and 6 have been obtained from a fit to the data. 
The three lowest momentum bins are most relevant to D*ii'. The excellent agreement between the 
corrected MC and the data can clearly be observed in this figure. 

The quadratic terms (S) may arise from differences in D* polarization; they are large (~ 0.5 
and 0.7) in bins 2 and 3 where we expect the D* spectrum to be dominated by highly polarized 
few-body B meson decays. The linear terms can only be attributed to differences in efficiency. 
The linear terms for the first three bins are small; i.e., they are 0.037 it 0.16, —0.023 it 0.024 and 
-0.016 ±0.009. 

5.2 Final state radiation 

The program PHOTOS |16j is used to model the effects of final state radiation (FSR). PHOTOS 
uses QED to second order in aem (up to two FSR photons can be produced) and is known to 
provide an accurate simulation of the FSR effects. 

5.3 Signal 

To simulate the signal we use Eq. for the distribution of the decay products. The MC samples 
are generated with the default parameters Ri = 1.180, R2 = 0.720 and = 0.920 P^ . 

5.4 Other semileptonic decays 

A major source of background is other semileptonic B decays. Aside from branching fractions 
for the decay modes Diu and the signal mode D*iu, only that for the mode B^Dilu has been 
measured [llj . Thus, the other branching fractions used in the MC simulation are based on theory. 
Until such time as we can measure them, we must accept quite large possible variations in their 
numerical input values in our simulation. 

6 Analysis Method 
6.1 Fitting 

Our basic approach to extracting the D*^v form factors is to perform an unbinned maximum 
likelihood fit to the full four-dimensional (4D) distribution (usually referred to as the probability 
density function or PDF in the context of likelihood fitting) specified by Eq. We parameterize 
the form factors in terms of the HQET inspired parameters Ri, R2 and as described in Sec. [2 
Since theoretical predictions differ on the w dependence of i?i and R2 (see Eas. ()17l23|) ), we first 
give as our baseline result the values obtained by treating them as parameters which are constants 
over the entire range of w {i.e. setting the coefficients of the {w — 1) and {w — 1)^ terms to zero and 
fitting only for the constant terms). This simplifies comparisons between experiments. However, 
we also show how the results vary when using the ii;-dependencies suggested by the full predictions 
of Eqs. JHESl). 

In addition to assuming the form of the theoretical PDF, we must account for the effects of 
resolution and efficiency on the kinematic variable distribution. As our Monte Carlo samples are 
limited in size and do not allow us to adequately map the resolution and efficiency as a function 
of -u;, cos 9i, cosOy and Xj we resort to the approximations derived below to allow us to correct for 
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the effect of efficiency and resolution on tlie fitted parameters using only integrals of this efficiency. 
We must also account for the residual background. 

To account for the efficiency and resolution effects we adapt the approach first employed in our 
angular analysis of the decay B J/tpK*^^ (for more details see also |14j). The PDF including 
resolution and efficiency {F) is given by 

F{x;fi) = J dxe{x)G{x;x) F{x;fi) (28) 

where x represents the true variables {w, cosOi, cosOy and x)^ ^ are the observed values of the 
variables and fi represents the parameters {Ri, R2 and p^) that determine the form factors. The 
efficiency £{x) is the fraction of events with parameters x that will be - on average - detected 
and G{x; x) is the probability that an event with true parameters x will be reconstructed with 
parameters x. 

The log of the likelihood L that we need to maximize is given by 

InL = ^In =Y.\nF{x,;fi) - N,,,,,^, x lniV(/.) (29) 

where the integral N{fj,) = J dxF{x; fj.) is needed to properly normalize the likelihood in the presence 
of imperfect acceptance. The sum is over our data event sample. A'^gignai is the number of signal 
events in the data sample. 

Now consider the following trivial modification (multiplication by 1 = F (x; fiuic) / F {x; fi^c)) , 
where fimc is the parameter set used to generate our Monte Carlo sample) of Eq. ^ 



F{x; ^i)= fdx e{x) G(i ; x) F{x; /i^c) x (30) 

If fi is not too different from then the ratio F (x ; fi)/F{x; fi^ac) will not vary much across 
the range where the resolution function is much different from zero. In this case we can use x to 
approximate x, which allows us to pull this ratio out of the integral to obtain 

F{x-f,) « -^^^ I dxe{x)G{x■x)F{x■^l^,) = F{x■,^,)^^^^. (31) 

F{x;flrnc)J F{x;^Xmc) 

When substituted into the expression for the log-likelihood this gives 

InL = ^\nF{xi;ii) - ^ lnF(xi; ^mc) +^\-n.F{xi; p,ri,c) - iVsignal lniV(/i). (32) 

Since terms that are independent of the fit parameters (constant terms) do not affect the point 
at which the maximum will be found, all the sums that depend on /_fmc can be dropped. The 
^-dependent piece has been factored from these constant terms. We are left with a likelihood 
that depends only on the theoretical PDF F and on the integral over the resolution and efficiency 
functions. Thus, there is no need for a detailed understanding of the efficiency and resolution as 
a function of the kinematic variables. All we need is a method for evaluating of the normalization 
integral N{^) 

We can use the technique of Monte Carlo integration to evaluate the integral N{fi). We have 
iV(,) . j dxFix; ,^.) X . ^ ^ ^^——^ (33) 
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where we have used the fact that the Monte Carlo generates events proportional to F{x]^^c) to 
approximate the integral by a sum over MC events divided by the number of events generated, 

We call this technique the 'resolution-efficiency correction' (REC) method. It can be shown to 
account for both the biases produced by resolution and efficiency if the Monte Carlo sample used 
was generated with parameters sufficiently near the final fitted values. As will be seen in the next 
section the small residual bias is also correctable. 

The last remaining step is to handle the background. Ordinarily, we would add a PDF B[x) 
that models the background. The PDF F would be replaced with 

fF{i;ix) + {l- f)B{i) (34) 

where / is the signal fraction. However, since we do not have a form for the background before 
acceptance, we can not achieve the factorization of the parameter dependence from the efficiency 
and resolution functions that leads to Eq. (|32p. To avoid this problem we subtract the background 
directly in our likelihood sum rather than adding it to our PDF. We replace our log-likelihood 
function with the following 'pseudo-likelihood' 

In^ = El^^(4L;/^) - E^Sgdlni^(x£gd;/^) - A^signailniV(/i)) (35) 

i i 

where the first sum is over the data and the second is over a Monte Carlo sample representive of 
the background. This method is not as statistically optimal as using a proper likelihood - however, 
it is still unbiased. 

The weights w^kgd ficcount for the normalization difference between the background in the data 
and in the Monte Carlo. They are computed in the manner indicated by Eqs. H41|) and ()42() in Sec. 
16.21 below. We call this technique for handling the background the 'direct unbinned background 
subtraction' (DUBS) method. 

The last term in Eq. (jSSJ) involves a sum (Eq. (jSSI)) over the Monte Carlo data sample. This 
could be quite computationally intensive. However, we are able to speed up this computation by 
transforming the sum over events into a sum over 'moments.' 

Since our PDF can be written in the following form, 

F(i;^i) = ^^„(/i) X ^„(x), (36) 

a 

i.e., as sum over a product of terms depending only on the fit parameters and terms depending 
only on the kinematic variables, we can define moments Mq by 

"^ = ^T.M^. (37) 

gen ,^ ^ \'^;H'mcJ 

where the sum is over reconstructed MC events, i.e., the same sum that defines N{fi) in Eq. ()33() . 
This allows us to write N{fi) as a sum over moments: 

N{fi) = Y,AMxM^. (38) 

a 

The moments can be computed once before fitting. In the fit, the time consuming event sum is 
replaced with the sum over moments. In our case taking the expansion of to order [w — 1)^ we 
have 18 moments to compute and sum. 
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6.2 Corrections and errors 



The approximations outlined above (the REC and DUBS methods) provide us with a fast and easily 
implemented fitting procedure that uses the MC sample without any need to extract a detailed 
model of the efficiency and resolution functions. However, these advantages do not come without a 
price: we must make corrections for the small residual bias and we must account for the uncertainty 
introduced by the REC and DUBS procedures. 

The REC method would be unbiased if the parameters fimc used to generate the Monte Carlo 
sample turned out to be the same as the true parameters. This is, of course, in practice almost never 
the case, so we must ex post facto introduce a small correction. To do this we need the derivatives 
of the fitted parameters with respect to the input parameters. We obtain these by reweighting the 
Monte Carlo sample to correspond to parameters that deviate by ±0.1 from the parameters used in 
generation {Ri = 1.18, i?2 = 0.72, = 0.92) and observing the differences produced in the fitted 
parameters. This yields the full 3x3 matrix of the derivatives d/ifitted/^^/^true- 

Given the derivative matrix we can relate the results of our fit to the corrected results (^corrected) 
as follows 

(^fitted - f^mc)= , ^ (/^corrected Mmc)' 

(39) 

To solve for the corrected parameters we need only invert the derivative matrix and multiply it 
against the LHS. The size of this correction is ~ 0.05 for Ri, ~ 0.07 for R2 and ~ —0.04 for p'^. 

Now since we use the theoretical PDF in our fit without explicit smearing to account for 
resolution, we have also underestimated the error. While the REC method corrects the values of 
the parameters for resolution effects (to a first approximation), it does not account for the loss 
of resolution due to our imperfect estimates of the four kinematic variables (which arise largely 
from not being able to obtain the exact B rest frame, due to the missing neutrino, but also from 
detector smearing). However, the error matrix -Efittcd for Ri, R2 and p'^ can be corrected using the 
derivative matrix, as follows 

i^fitted ^ ( ^) " i^fitted ( ( f^^) " V . (40) 

This increases the error estimates by ~ 12%. 

We call this the 'bias map correction' (BMC) procedure, which refers to both the modifications 
of the central values, and to increase in the estimated statistical uncertainty. 

In addition to the BMC error increase, there are contributions to the error that are not ac- 
counted for in the fit. Two of these are Monte Carlo statistical in nature: the error from the 
DUBS subtraction and the error from the Monte Carlo integration used to evaluate N{p) in the 
REC procedure. There is also a contribution to the statistical error from the fluctuations in the 
background that is not accounted for in the fit. The first two can be reduced by increasing the size 
of the Monte Carlo sample, but the latter error is irreducible and must be included regardless of 
the size of the Monte Carlo sample. 

We have developed procedures for evaluating these three errors using appropriate sums of the 
derivatives of the likelihood that can be evaluated using the Monte Carlo sample. The formulas 
used are collected in Appendix A. 



22 



6.3 Background 

We rely on Monte Carlo to model the background. However, in order to subtract the background 
it must be properly normalized. To do this we adopt the background estimates from our Vcb 
analysis [Sj. The Vcb analysis estimates the background by fitting the cos and 6m distributions. 
The Vcb analysis fits the cos 9 by distributions for seven different background contributions (obtained 
either from MC or data sidebands, see jS| for a description of the seven categories). For our 
purposes we lump the background into two categories - peaking (candidates that contain a properly 
reconstructed D* and hence peak in the 5m distribution) and combinatoric (candidates in which 
the D* is a fake made from an incorrect combination of particles). Combinatorial is one of the 
seven categories above; we obtain the peaking fraction by summing over the other six. 

The fractions of peaking and combinatorial background are determined in ten it;-bins. Thus, 
we do not rely on the MC to model the w-dependence of the background, but take it from data. 

For each w-hin we compute weights to apply to each event from the MC background sample 
before subtracting it. For the peaking component the weight for event i is 

^"peaking = /peaking (w^) X A^data X (41) 

where the Wi are weights we apply to correct for difference in the numbers of B^B^, B~^B~ and 
continuum Monte Carlo events available and A'^jata is the number of events in our signal window. 
The sum is over peaking background events {k) that are in the w-hin where the w of the candidate 
falls. 

This procedure guarantees that the number of peaking background events subtracted will be 
equivalent to the number (/peaking (i*^) x A^data) estimated to be inside the signal window. 
A similar prescription applies for the combinatoric background with weights defined by 

^cLh = fcomhH X A^data X (42) 

where the events k are now selected from the combinatoric category. 

7 Results 

For the baseline result we perform the fit taking Ri and R2 to be independent of w. All results are 
preliminary. 
We find 

Ri = 1.328 ± 0.055 ± 0.025 

R2 = 0.920 ± 0.044 ± 0.020 (43) 
p2 = 0.769 ± 0.039 ± 0.019 

where the first error is statistical and the second Monte Carlo statistical. The BMC fEans. H39ll40() ^ 
has been applied to both the central values and the the statistical error (it shifts the former and 
expands the latter). Systematic uncertainties are discussed in Sec. El below. 

The errors are highly correlated. The error matrix for the full statistical error for Ri, R2, and 

(including Monte Carlo) is: 
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0.003664 
-0.001696 
0.001213 



-0.001696 
0.002342 
-0.001634 



0.001213 
-0.001634 
0.001850 



The correlations are 

PR1-R2 = -0.58 

PR,^p2 = +0.47 (44) 
PR,.p2 = -0.79. 

As we do not yet have enough sensitivity to fit for u;-dependence of Ri(w) and R2{w), we 
consider instead the effect of the theoretically predicted dependence on the result. Parameterizing 
this dependence as follows 

Ri{w) = Ri + ai{w - 1) + Pi{w - if (45) 

R2{w) = R2 + a2{w - 1) + I32{w - if (46) 

and inserting these t(;-dependent forms into the PDF (with fixed ai,f3i from the theoretical predic- 
tions) and fitting for the constant terms Ri and R2 we find the results given in Table ^ 



Table 1: Dependence of form factor parameters on theoretical assumptions about 
slope (a) and curvature (/3) of Ri and R2 w-dependence. See Eqns. and 



Reference 


ai Pi a2 132 Ri R2 


NeubertP 


-0.22 0.09 0.15 -0.04 1.37 0.88 0.75 


Caprini-Lellouch-Neubert 


-0.12 0.05 0.11 -0.06 1.35 0.89 0.76 


Ligeti-GrinsteinjS] 


-0.10 0.0 0.09 0.0 1.35 0.89 0.75 


Baseline 


0.0 0.0 0.0 0.0 1.33 0.92 0.77 



Taking into account these theoretical variations in Rii^w) and R2{'w) yields slightly larger values 
for Ri and slightly smaller values for R2. The prediction by Neubert produces more deviation than 
the more recent calculations. 

8 Goodness-of-fit 

As these results are obtained from a maximum likelihood fit, we need some method of assessing 
whether the results of the fit actually reproduce the distribution of the kinematic variables in 
the data. Since we do not have an explicit form for the acceptance-corrected PDF (F) of the 
four reconstructed variables w, cos 61, cos By and X; we reweight the MC sample to construct the 
distributions expected from our measured parameters. That is, the contribution of the signal to a 
bin is 

r^signal = ^lignal (47) 
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where 

^iignal = (1 - /peaking (W^) " fsignsd{w)) X iVdata X (48) 

and in this case Wj = J^^^'^^ is the weight needed to modify the distributions from those generated 

* F{x;llrnc) to J t) 

with fimc {Ri = 1-18, R2 = 0.72, and = 0.92), to those we obtain from this analysis. 

For the background we use the same weighting procedure (see Eqs. (|1T|) and (|42|) ) used in the 
fit. Using these weighting procedures the normahzations of the data and reweighted distributions 
match by construction. 

We consider two types of goodness of fit: (a) five one-dimensional distributions (the projections 
w, cosOi, cos9v, and x plus the distribution of CM lepton momentum p|) and (b) a binned 
based on 6x6x6x6 (a total 1296) bins. For a comparison we also give the unweighted results 
which corresponds closely to those using the CLEO parameters 0- 

In Figs.[3nia comparison of the projections in the kinematic variables and between the Monte 
Carlo and the data can be seen. The plots in Fig. [7| show the result for the default parameters, 
while those in Fig. |H1 give that obtained by reweighting by our parameters, and Fig. IHl shows both 
for the p£ spectrum. The fit is to a constant from which we read off the x^- 

Close inspection will reveal that the ratio plot projections in Fig. |H1 show better agreement to 
the line at unity than those in Fig. Q but the improvement is more clearly seen numerically in 
Table [21 which gives the x'^ for fitting each plot to a constant. The constant is always consistent 
with unity as it must be. In every case the agreement between MC and data improves when we 
use our result - sometimes substantially. 

Table 2: x^ and x^-pi'obability for kinematic variable projections and lepton momen- 
tum; the number of bins in these histograms is 10. The number of degrees-of- freedom 
is 9, since we have forced the normalization. 



variable 


X^ (prob.) default 


X^ (prob.) BaBar 


w 


14.08 (0.120) 


13.69 (0.134) 


cos 9i 


17.12 (0.047) 


7.626 (0.572) 


cos By 


40.93 (0.000) 


17.81 (0.0374) 


X 


8.133 (0.521) 


7.082 (0.623) 


* 

Pi 


17.88 (0.037) 


7.316 (0.604) 



The value of the four-dimensional binned goes from 1274.04 with default MC parameters 
to 1232.23 with our parameters - an improvement of ~42 units of x^ . The x^ per bin is slightly 
smaller than unity. It is not really proper to interpret this as a probability as there are about 200 
empty bins which contribute nothing to x^- Nonetheless it is clear that the reweighted MC follows 
the distribution of the data quite well and that it models this decay (the single largest B branching 
fraction) much better than the default MC. 
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Figure 7: Data (points) overlayed on Monte Carlo (histograms) for all four kinematic variable 
distributions for default values of the parameters The plots below the overlays are the ratios of the 
data to MC distributions. 
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Figure 8: Data (points) overlayed on Monte Carlo (histogram) for all four kinematic variable 
distributions use our measurement of the parameters The plots below the overlays are the ratio of 
the data to the MC distributions. 
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9 Systematic Studies 



The systematic uncertainties on the three principal parameters are summarized in Table The 
dominant systematic errors arise from the MC simulation: that is, from how well we understand and 
simulate the detector performance in terms of resolution and efficiencies, in particular the efficiency 
for the reconstruction of low-momentum charged pions from D*+ decays. Further, how well we 
model the signal and background event generation, e.g. how close the branching fractions in the 
event generator are to those of the real world, affects the distribution of background we subtract. 
The w-dependence of the background, however, is not taken from the MC but is measured in the 
data and the uncertainty in this measurement also contributes to the systematic errors. 

9.1 Detector Simulation 

Extensive studies of the simulation of the detector response, including careful examination of track 
reconstruction efficiencies and particle identification, have been performed using selected data con- 
trol samples. Adjustments for known simulation deficiencies are used in investigating and evaluating 
the systematic errors. Form factor measurements are insensitive to overall normalization errors. 
Thus differences in the efficiencies that are independent of the fit variables do not affect the results, 
but variations of differences as a function of these variables are of concern. 

To assess the uncertainties due to differences in shape rather than normalization we vary the 
dependence on the efficiency corrections, reprocess the MC samples, redo the fit to the data, and 
take the difference to the results obtained with the nominal MC simulation as an estimate of the 
systematic error on the parameters. This procedure is repeated for every source of systematic 
uncertainty. The individual uncertainties are added in quadrature to obtain the total systematic 
error. 

9.1.1 Charged Particle Tracking 

The difference between the tracking efficiency for electrons, and charged kaons and pions from D 
decays decreases roughly linearly as a function of momentum. We vary this linear dependence 
on the particle momentum in the MC simulation and observe no significant change in the results. 
Consequently, we assign no error from this source. 

9.1.2 Slow Pion Reconstruction 

The efficiency for reconstructing the low-momentum charged pion (tTs) from D*^ is a major source 
of systematic errors. Because of the small energy release in D* decays, this pion is emitted in 
the same direction as the parent D* and its momentum is less than 400 MeV/c in the laboratory 
frame. Since w = E£,*/mD* in the B rest frame, the tt^ momentum is correlated with w and thus 
its momentum dependent efficiency impacts the measurement of p^. 

The uncertainty due to the low-momentum tracking efficiency is evaluated differently from other 
tracking errors because such low-momentum tracks do not traverse the whole drift chamber. Their 
detection and measurement depends mostly on the silicon vertex tracker. To study this efficiency 
as a function of p-,^^ we use a large set of D*"^ — > D^irf decays selected from hadronic events and 
measure the distributions of the helicity angle of vr^ (^y) in the D*~^ rest frame as a function of 
the D*~^ momentum. Fig. 1^1 in Sec. 14.21 compares the cosOy obtained for data and MC simulation. 
The small size of the linear contributions to the correction factors (Eqn. I27|) is encouraging, but 
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there is enough room in the error that we still need to include the vr^ efficiency in the systematic 
error. 

We parameterize the effective vr^ efficiency as function of its momentum using the form 



with po being the threshold momentum and /3 controling the rapidity with which the efficiency rises 
above threshold. For p^^^ < pQ we set e to zero. We fit the data and MC helicity angle distributions 
for (3, po and the coefficient of the physically allowed cos^ Oy term to obtain efficiency functions for 
the data and the MC. The helicity method only determines the relative momentum dependence of 
the efficiency, so the normalization Emax must be determined separately. Since, normalization does 
not matter for this analysis we simply set it to unity. 

To assess the systematic uncertainty due to the vr^ efficiency we weight the MC simulation by 
the ratio of data to MC functions and assign the observed shifts in the fitted values for Ri , i?2 and 

as systematic errors. Not unexpectedly, is most sensitive to this efficiency since it describes 
the shape of the t(;-distribution. 

9.1.3 Charged Particle Identification (PID) 

Using data and MC simulated control samples we have tabulated the difference in particle identifi- 
cation efficiency for data and MC. These tables provide the correction factors in bins of momentum 
and angle. Since this analysis is most sensitive to efficiency variations with momentum we average 
the tables over the angles to obtain corrections as a function of momentum. For electrons the 
correction factors vary from 0.991 to 1.008 over the momentum range from 1.2GeV/c to 2.5GeV/c. 
We assess the impact of the uncertainty in these corrections by approximating their momentum 
dependence by linear functions and vary the sign of the small slope of these functions. The ob- 
served deviations from the default fit are ARi = 0.0064, AR2 = 0.0052 and Ap"^ = -0.0016 for the 
positive slope and -0.0032, -0.0031, -1-0.0009 for the negative slope. We take half of the difference 
as the systematic error from this source. Since the momentum dependence is not a monotonic 
function, this procedure overestimates the uncertainty. 

For kaon identification we employ the same procedure. The observed variations are significantly 
smaller. 

The probability of misidentifying charged hadrons, vr^, K^, p^, as electrons is very small, less 
than 0.2% in the momentum range 1.2 — 2.5GeV/c. Since a variation of the peaking background 
by 9% results in a very modest change in the fit results, and since the fraction of this background 
originating from hadrons misidentified as electrons is very small, we conclude that the uncertainty 
in the hadron misidentification rate should be negligible. 

The misidentification rate of pions as kaons ranges from a few tenths of a percent to almost 5%. 
However, pion misidentification is well simulated by the MC and thus should have little impact 
on the fit results. Furthermore, the main consequence of pion misidentification is to increase the 
combinatorial background. Since we estimate the combinatorial background from a fit to the 
measured 5m distribution, we are not dependent on the MC to assess the size of this background. 
We conclude that the uncertainty in the pion misidentification rate has little impact on the fit 
results. 




(49) 
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9.2 Event Simulation 



9.2.1 Final State Radiation (FSR) 

Final state radiation, primarily from electrons, lowers the momenta and to a lesser degree changes 
the angles of detected particles. Though a physics effect, FSR acts much like a resolution - it 
smears the kinematic variables. We simulate the emission spectrum of radiative photons using 
PHOTOS ^ni; so the REG method corrects for it, to the extent that PHOTOS models it correctly. 

To test the sensitivity to FSR we evaluate the shifts in the fitted values of Ri, R2 and p'^ 
between fits done with and without FSR corrections. We assume an uncertainty of 30% in the 
simulated photon emission and thus take ~ 1/3 of the observed shifts of 0.0129, 0.067 and 0.0039 
as an estimate of the systematic uncertainty. 

9.2.2 Background Simulation 

We divide the background sources into two categories: (1) peaking background for which the D* de- 
cay has been correctly reconstructed and which contributes to the peak in the 5m = nix-K+w- —^Ktt 
distribution and (2) combinatoric background for which the D* has not been properly reconstructed 
and thus does not contribute to the peak in the 6m distribution. 

Background mixture 

Our modeling of the peaking and combinatorial background depends on the not very well known 
branching fractions for the mixture of semi-leptonic B decay modes that make up the background. 
To estimate the uncertainty associated with these branching fractions, we vary their values and 
observe shifts in the fit parameters compared to the nominal values. In this process we keep the 
total background fractions as determined by the cos 6by ht (see Sec. 16. 2|) unchanged. We vary 
most modes by 60%. For the measured mode B^Dilv we vary only by the ~ 30% measurement 
error JJ. In the case of D*lu there are contributions from badly reconstructed signal events. We 
assume a 20% uncertainty in the signal branching fraction to account for the large differences in 
the currently available measurements. 

In TableOlwe list these branching ratio variations as well as the effect of varying the contribution 
from e"'"e~ — > cc events. We take half of the observed variation in the fit parameters as an estimate 
for the systematic error. 

The total error for the three fit parameters is obtained by adding the errors due to each contri- 
bution in quadrature. Ri is the most sensitive to the mixture of decay modes of the background 
subtraction. 

Dependence of the Background on w 

The tt;-dependence of the background estimate is taken from cos^^y fits performed for each bin 
in [H]. We fit the it;-dependence of the peaking and the combinatorial backgrounds to a second 
order polynomial centered at the middle bin. We use these polynomials, fpcak\ng{w) and fcomhiw), 
to compute the weights (see Eqs. (|1T|) and H42() ) that normalize our background subtraction. 

To estimate the error from the w-dependence we vary the slopes of /peaking (^^) and fcomb{w) 
polynomials by itlcj and refit for R2 and with the altered background normalizations. We 
take half of the total variation from each background type as its contribution to the systematic 
error. 
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Total peaking and combinatoric background fractions 

Another source of error is the uncertainty in the normalization of the peaking and combinatoric 
background fractions. To estimate this error we vary the fraction /peaking by ±12% and /comb by 
±15%. We have increased these uncertainties beyond the statistical uncertainites established by our 
Vcb analysis |S] to account for slight differences between the way the analyses define backgrounds. 
For /peaking this yields errors of 0.0186, 0.0075, 0.0014 for the three parameters, indicating that this 
is a significant error source for Ri. For /comb we find error estimates of 0.0063, 0.0034, 0.0080. 

MC/data side-band comparison 

The distributions of the kinematic variables for MC and data agree well in the 6m sideband 
region used to estimate the combinatoric background. But since the combinatoric background 
comprises about a third of the total background under the peak, the small differences in the shapes 
of the distributions could introduce an error in the background subtraction process (when the MC 
is used to subtract the residual combinatoric background from the data, see Sec. 16. If) . To estimate 
the impact of this effect, we first take the ratios of data to MC in the sideband region and then fit 
the distributions with polynomials with the result shown in Fig. 1101 

We then use these functions one at a time to multiply the combinatoric background from MC 
before it is subtracted from the data to prepare the sample for fitting. We then carry through the 
fits and then take the differences in the form factors we obtain from these with the form factors 
we obtain from the fits with the unaltered background. This procedure yields the results shown in 
Table H 

The differences we find are small. The largest is from using the function for w, from which we 
find Ap^ 0.006. Since in the end we take the tf-dependence of the background from the data 
and the deviations due to weighting the angular distributions are very small, we add nothing to 
the systematic uncertainty from this check. 



9.3 Summary of Systematic Errors 

The systematic errors are summarized in Tablel^l While the total error remains statistics-dominated, 
we are approaching the systematic limit. The addition of other decay modes with higher mul- 
tiplicity and higher backgrounds is likely to increase our sensitivity to the various background 
sources. Enhancing our understanding of these backgrounds will be critical to improving the er- 
rors. Measurements of the higher mass semileptonic decay modes are sorely needed. 

10 Conclusions 

We have measured the form factors Ai, V and A2 in terms of the HQET-inspired parameters Ri, 
R2 and p^. Note that all results are preliminary. 
The baseline result including systematic errors is 

i2i = 1.328 ±0.055 ±0.025 ±0.025, (50) 
R2 = 0.920 ± 0.044 ± 0.020 ± 0.013, 
p2 = 0.769 ± 0.039 ± 0.019 ± 0.032 
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where the first error is statistical, the second MC statistical, and the third systematic. This baseline 
result is obtained neglecting any possible ^-dependence of and R2{w). It agrees well, but 

not perfectly, with theory at w = 1 (see eq. H17() or (I2U() ). We achieve a considerable improvement 
over the CLEO result of i?i = 1.18±0.30±0.12, R2 = 0.71±0.22±0.07 and = 0.91±0.15±0.06 

We do not yet have the sensitivity to independently establish the w-dependence of Ri{w) and 
R2{w), but we can use theoretical estimates to make comparisons. If we compare the predictions of 
Caprini, Lellouch and Neubert |1] for Ri and R2 to the result we obtain using their w-dependence 
(see Tabled, we find ARi = 1.36 - 1.27 = 0.09 (1.6cr statistical) and AR2 = 0.87 - 0.80 = 0.07 
(1.4(7 statistical). Given the undoubted presence of theoretical errors, this is reasonable agreement. 

This measurement allows us to make a considerable improvement in the size of the current error 
in Vcb measurements. For our measurement of Vcb the error contributed by uncertainty in Ri and 
i?2 drops by a factor ~ 4 from that obtained using the previous measurement from CLEO [2] . The 
effect on the total uncertainty is to reduce it by ~ 30%. 

A considerable improvement can also be obtained in measurements of the lepton end point 
spectrum. The systematic error on the branching fraction for decays with a lepton in the momentum 
range 2.1 — 2.6GeV/c is reduced from 5% to 1.7% which corresponds to an improvement of ~ 8% 
in the total statistical plus systematic error. 

In addition we have demonstrated useful approximations to the maximum likelihood method 
that allow us to cope with the limited size of the Monte Carlo samples available to modern high 
luminosity experiments. We have also developed the procedures needed to evaluate the corrections 
and additional uncertainties due to these approximations. These methods are not unique to 0*1^ or 
J/ipK* , but could be applied to any analysis that needs to cope with a complex multi-dimensional 
acceptance and resolution problem. 
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0.8 U 



Figure 9: Data/Monte Carlo comparison for p|-distribution (data are points, histograms are MC). 
MC used for the top plot uses default parameters and MC used for the bottom is reweighted to 
our parameters. 
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Table 3: Systematic errors due to the uncertainties in the branching fractions of 
semileptonic B decay modes. The branching fractions hsted are the default values 
used in the MC simulation. The isospin related modes of neutral and charged B 
mesons are varied together by the factors indicated. 



1 lor'PiT' ^/li^riti 


IVlV^ Ul dlidillie, iidjULUJli I /U 1 


\/QT'lQ'l'li~m / /r\ \ 

V di IcLLlUli \ /C / 








B^D*lv 


5.6 


20 


0.00052 


0.00044 


0.00027 


B^Diu 


2.10 


20 


0.0013 


0.00037 


0.0002 


B^Doliy 


0.20 


60 


0.00020 


0.00026 


0.00010 


B^Diiu 


0.56 


30 


0.0087 


0.0.0024 


0.0016 


B^D[£u 


0.37 


60 


0.012 


0.0062 


0.0044 


B^D2£v 


0.37 


60 


0.00095 


0.0025 


0.0017 


B^D*7t£u 


0.30 


60 


0.0036 


0.00087 


0.00071 


B^DttIv 


0.9 


60 


0.0022 


0.00092 


0.00061 


e+e~ cc 


NA 


20 


0.0011 


0.00034 


0.00040 


Total 






0.016 


0.0073 


0.0051 
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Figure 10: Ratio of data to MC for the four kinematic variables in the high sideband region data, 
and the best polynomial fits. 
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Table 4: Changes in the fitted parameters for reweighting of the MC combinatoric 
background distributions in the four kinematic variables, as shown in Fig. IIUI 



Reweighted distributions 


Ri 


R2 




w distribution 


-0.002 


0.0 


0.006 


cos^^ distribution 


0.001 


-0.002 


-0.001 


cos^v distribution 


0.002 


-0.003 


0.001 


X distribution 


0.004 


-0.002 


0.001 



Table 5: Summary of the estimated Systematic errors. Negligibly small contribu- 
tions have been omitted here, they are discussed in the text: tracking, FSR, and 
5m MC/data differences. 



Error source 








Lepton-hadron track efficiency 


0.005 


0.004 


0.002 


Slow pion track efficiency 


0.003 


0.0002 


0.011 


PID misID (lepton, kaon) 


0.005 


0.004 


0.002 


Peaking background normalization 


0.019 


0.008 


0.001 


Combinatorial background normalization 


0.006 


0.003 


0.008 


Background composition (branching fractions) 


0.016 


0.0079 


0.0051 


tf-dependence of background 


0.001 


0.001 


0.028 


Final state radiation 


0.0043 


0.0023 


0.0013 


Total Systematic 


0.027 


0.014 


0.033 
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Appendix A: Evaluation of additional errors 

Here we gather the equations used to evaluate the errors not included in our likelihood fit. 
Specifically, we provide equations for the errors due to the Monte Carlo integration of the REC 
method, that due to the background subtraction of the DUBS method, and that due to the fluctu- 
ations of the actual background initially left unaccounted for when the DUBS method is used. 

Full details and derivations can be found in the thesis of M.S. Gill[14j. 

REC Monte Carlo integration error 

The REC Monte Carlo integration error is given by: 

Erec = EfitDE].^ (51) 

where Efn matrix is the error reported by the fitter. D \s a, matrix of derivative sums. 
D is composed of three pieces. To lay them out we need to define some terms. 
The sums of the weights Wi = F{xi; fi)/F{xi; fj-mc) used in the integral: 

W{i,) = Y.Wi (52) 

where sums over MC events used to evaluate the normalization integral N^n). 
The vector of the derivatives weights: 

^ dm _ ( dwj dwj dwi \ 

where Fi = F{xi; p). The derivatives are computed numerically though an analytic compuation is 
possible. 



The sum of weights squared: 
The W derivative vector: 



which also allows us to compute ^^^^ = w%:- 
Weight X derivative vector: 

WdW = ^ X Ai (56) 

Derivative x Derivative matrix: 

dWdW = (^iM) (57) 



The three components are give by: 

Di = r]X dWdW (58) 

((^) ------ (^)l (-) 

where rj = Ngignai/W^ is a normalization factor. 

D is just the sum of these three, i.e., D = Di + D2 + D3. 
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DUBS error 

The DUBS error can be computed from sums over the Monte Carlo sample used in the back- 
ground subtraction. In this case the weights (wi) are those used to weight each type of background 
to obtain the correct normalization as described in Sec. 16.21 fsee eqs. H41() ()42() ). 

In this we need the vector of the derivatives of the PDF from which we compute the 

'sensitivity' matrix S: 

S^u.s = Y.^.[-^)[-^) (61) 
Given S the DUBS error is just: 

Edubs = EfitSdubsE^fit (62) 



Background error 

We also use the DUBS sample to estimate error from background in our signal sample. The 
result is analogous to eq. (|HT|l accept that wf is replaced by Wi. That is, we have: 

^..^E».(^)(^)' (as) 



and 



Eback — EfitSbackEji^ (64) 
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